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The relativistic accretion rate of dark matter by a black hole is revisited. Under the assumption 
that the phase space density indicator, Q — poc/o^,, remains constant during the inflow, the derived 
accretion rate can be higher up to five orders of magnitude than the classical accretion formula, valid 
for non-relativistic and non-interacting particles, when typical dark halo conditions are considered. 
For these typical conditions, the critical point of the flow is located at distances of about 30-150 
times the horizon radius. Application of our results to black hole seeds hosted by halos issued from 
cosmological simulations indicate that dark matter contributes to no more than ~ 10% of the total 
accreted mass, confirming that the bolometric quasar luminosity is related to the baryonic accretion 
history of the black hole. 
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I. INTRODUCTION 

The presence of supermassive black holes (SMBHs) in 
the center either of elliptical galaxies or bulges of spiral 
galaxies seems presently to be a well established fact [l[ . 
However, in spite of a large number of studies performed 
in the past years, the origin and evolution of these objects 
remain quite uncertain. Since SMBHs cannot be formed 
as a consequence of the evolution of stars having compa- 
rable mass, they are generally supposed to evolve from 
"seeds" , which grow by accretion inside the host galaxy. 
Moreover, the observation of quasars already at z ~ 6.5, 
when the universe was only ~ 0.85 Gyr old, implies that 
the process of growth (accretion) must have been quite 
efficient at those early times. 

Seeds could be intermediate mass (10 3 -10 4 Mq) black 
holes formed in the collapse of primordial gas clouds 
@) S 0] or in the core collapse of relativistic star clusters 
formed in starbursts, which may have occurred in the 
early evolution of galaxies (5j. Another appealing pos- 
sibility is the formation of black hole seeds with masses 
in the range 100-300 Mq, resulting from the evolution 
of primordial massive stars [|| . These stars would have 
been formed at z ~ 15-20 in the high density peaks of the 
primordial fluctuation spectrum and their high masses 
would be a consequence of a very inefficient gas cooling 
at zero metallicity [7|. Besides forming black hole seeds, 
these massive stars could also have been the source of UV 
photons, responsible for the reionization of the universe. 

The analysis of the bolometric quasar luminosity func- 
tion at different redshifts indicates that most of the ac- 
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creted mass by black holes should be baryonic in origin 
[li H OS- Although expected to be small, the contri- 
bution of dark matter to the black hole growth is still 
uncertain. Unlike baryons, the infall of dark matter into 
black holes does not produce electromagnetic radiation 
and thus, does not contribute to the quasar luminos- 
ity, although some energetic photons could be produced 
either by annihilation or decay, if dark matter is con- 
stituted by neutralinos or sterile neutrinos respectively. 
Moreover, dark matter particles constitute a collisionless 
fluid and, consequently, the accretion process should be 
less efficient than that expected for a dissipative fluid. 
Even so, different scenarios were envisaged in which dark 
matter could give a significant contribution to the mass 
accretion budget. If halos embedding galaxies are con- 
stituted by self- interacting dark matter [ll|, then black 
holes with masses > 10 6 Mq could be formed directly as 
a consequence of relativistic core collapse of halos, follow- 
ing the gravothermal catastrophe [12| . In this scenario, 
no baryons and no prior star formation are required to 
form the massive black holes, but the high interaction 
cross-section necessary for the mechanism to work re- 
mains a difficulty A different approach was considered 
in [l3| , where the dark matter accretion process was in- 
vestigated in a scenario in which stars are already formed 
in the core of the galaxy. Dark matter particles scatter 
on stars in the vicinity of the central black hole, produc- 
ing an important inflow as a consequence of diffusion in 
phase space and a mass growth scaling as Mbh oc i 9 / 16 . 
A similar mechanism was examined in [Tij but, in this 
case, dark matter particles are scattered by molecular 
clouds. According to the authors, this process can boost 
the growth of the mass of seeds from about ~ 5 Mq 
up to ~ 10 4 Mq, which then grow by Eddington lim- 
ited baryonic accretion. The growth of seeds under adi- 
abatic conditions was considered by (la |. who obtained 
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a mass-age relation only in the isothermal case but with 
the resulting mass-velocity dispersion relation differing 
from observations. However, they have obtained a better 
agreement with data by assuming that black holes grow 
proportionally to the host halo mass. 

In the present work, a new investigation on the dark 
matter accretion process into black holes is reported. Our 
approach gives an accretion rate significantly higher than 
that based on the well-known cross-section for capture by 
a black hole of non-relativistic and non-interacting par- 
ticles. We have performed some numerical applications 
of our formula to the growth of seeds inside dark mat- 
ter halos issued from cosmological simulations. These 
numerical calculations indicate that, if there is a coeval 
growth of seeds and their host halos, then dark matter 
contributes to no more than ~ 10% of the total accreted 
matter, giving a further support to the usual assumption 
that the bolometric quasar luminosity function is related 
to the baryonic accretion history of SMBHs. This pa- 
per is organized as follows: in Section II, the equations 
describing the accretion process are presented; in Sec- 
tion III, numerical examples are discussed and finally, in 
Section IV, our main conclusions are given. 



II. THE ACCRETION MODEL 

The relativistic steady spherical flow of a perfect gas 
into a black hole was first investigated by [16[ and later 
by [13) EM O H3] among others. Radiation transfer ef- 
fects were considered in the pioneering work by Thorne 
and collaborators ([Hl[13|) and more recently in [23L[2^|. 
A comprehensive review on the subject can be found in 
[25l | . It is worth mentioning that if the fluid has a non- 
zero angular momentum with respect to the black hole, 
then the spherical symmetry is broken. In this case, the 
accretion process is more complex since mechanisms of 
angular momentum transfer must be considered. Studies 
by [1^1 based on a "shell" accretion model indicate that 
the black hole growth factor is reduced if the fluid has a 
non-zero angular momentum. As we shall see later, our 
cosmological simulations do not have enough resolution 
to describe adequately the central region including the 
black hole influence radius. Thus, in the present investi- 
gation, we restrict our analysis to the spherical case. 

Past investigations have shown that for an adiabatic 
flow of a gas with 7 < 5/3, the only critical point of the 
flow which lies outside the horizon is the one correspond- 
ing to the Bondi solution. As we shall see below, a similar 
situation occurs in the case of a collisionless fluid. In our 
scenario, seeds appear around z~12 — 15 as a result of 
the core-collapse of massive stars @ and are located in 
the center of dark matter halos. Near the central black 
hole of mass Mth, we assume that the spacetime is de- 
scribed by the Schwarzschild metric and that a steady 
spherical inflow is set up inside the influence radius of 
the black hole, defined by r* = 2GM{,/ l /er 2 , where is 
the velocity dispersion of dark matter particles of mass 



m in the core of the halo. 

Under these conditions, the equations of motion are: 
the conservation of the mass flux 
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and the energy-momentum flux conservation 



(1) 



(2) 



where J k = mnu k is the mass-current density and the 
energy-momentum tensor is that of an ideal fluid, i.e., 
Tf = (P + e)u k m — PS k , with P, e, n, being respec- 
tively the proper pressure, the proper energy density 
and the proper number density of the dark matter fluid. 
u l = dx l / ds is the 4- vector velocity, whose non-null com- 
ponents are u = u 1 — dr/ds and u° = dt/ds. As in [l6l |. 
one obtains from these equations 
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where r g = IGM^jc 2 is the gravitational radius and k\ 
is a constant. Deriving eq. [3] with respect to the radial 
coordinate r and using the mass conservation equation 
(cq. [T]), one obtains a "wind" equation, i.e., 



dig u 
dlgr^ 

where 



and 



V 2 F(r, it)] + - 2V 2 F(r, u) 
2r 



F(r, u) 



(1-* 
r 



V) 



(4) 



(5) 



V 2 



dig (P + e) 
dig n 



1 



(6) 



The critical point of the flow occurs where both brack- 
eted factors in eq. [4] vanish simultaneously, namely, 
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After decoupling, dark matter particles constitute a 
weakly interacting fluid. In the matter era, at higher red- 
shifts, the universe evolves closely to the Einstein-de Sit- 
ter model. In this case, the distribution function f(u,t) 
of dark matter particles is a function of integrals of the 
characteristics of the Vlasov equation, i.e., 



f(u,t) = f(a(t)u) 
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where a(t) oc i 2 / 3 is the scale factor and u = (v — Hr) 
is the peculiar velocity. The distribution function must 
also satisfy the normalization condition 
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where the small contribution of baryons to the total mat- 
ter density was neglected. Under these conditions, the 
velocity dispersion of the peculiar velocity satisfies 
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and the quantity Q = p(t)/ < a 2 > 3 / 2 , an indicator of 
the phase space density, is conserved during the expan- 
sion of the universe. 

High resolution simulations of galaxy-size dark matter 
halos by (2|| indicate an increase of Q toward the center. 
Similar results were obtained by [30], who have derived 
a power-law variation, e.g., Q oc for cluster-size ha- 
los, with (3 quite close to the value found by [2{|, namely 
(3 ~ 1.87. No adequate explanation exists presently for 
such a power-law profile, although [3l[ suggests that it 
could be a consequence of an isotropic distribution func- 
tion produced in the halo core by processes of relaxation 
related to the radial orbit instability. In fact, [Uli^l as ~ 
suming that the evolution of halos is self-similar, derived 
the conditions required for the density profile as well as 
the phase density indicator Q be described by power laws. 
A density profile with a slope dlgp/dlgr = —2, typical 
of self-similar radial infall, can be consistent with phase 
space profiles either of the form Q — p/a 3 or Q = pj a 2 . 
Here we adopt the usual definition, which is supported 
by the results of different cosmological simulations. 

As structures begin to form as well as to accrete mat- 
ter and merger together, rapid variations in the gravita- 
tional potential contribute to the relaxation of the system 
through "violent relaxation" ([27[). This process leads 
to a smooth mass-independent distribution function as 
a result of the gravitational scattering of particles. Vio- 
lent relaxation produces a more mixed system in phase 
space, reducing the value of the coarse-grained distribu- 
tion function or, in other words, reducing the value of 
Q ((Hi)- Therefore, as a consequence of such a mech- 
anism, the velocity dispersion of dark matter particles, 
at late epochs, inside halos is not negligible as one would 
expect taking into account only the expansion of the uni- 
verse (eq. [TTj) . In this case, including the energy of the 
random peculiar motion, the total energy density can be 
written as 
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and, consequently 
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Notice that, as mentioned previously, the phase space 
density increases toward the center as Q oc r -1,87 , but 
one expects a convergence to a finite value close to the 
black hole influence radius. Once the inflow develops, we 
hypothesize that the quantity Q remains constant. This 
assumption is reasonable if no important phase space 
mixing occurs during the inflow. According to cosmo- 
logical simulations by [H, [H| , the phase space indicator 



Q estimated in the core of halos evolves in two distinct 
phases. An early and fast phase in which the decrease 
of Q is associated to virialization just after the first shell 
crossing, following the gravitational collapse of the halo. 
In a second late and long phase, a slow decrease of Q 
is observed, consequence of continuous matter accretion 
and merger episodes. Since the infall time scale of dark 
matter particles onto the central black hole is small in 
comparison with the characteristic time scale in which 
the phase space indicator varies, namely, Q/ \ (dQ/dt) |~ 
0.7 Gyr (typical value derived from the aforementioned 
simulations), we can suppose that the flow is steady but 
adjusted to the slowly varying conditions of the environ- 
ment. In this case, eq. [13] can be recast as 
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Notice that, in the above equation, we have defined 
Qn = Q/fn = n/a 3 and that the assumption of keeping 
Q constant during the inflow is similar to that of con- 
sidering the compression of a perfect gas with 7 = 5/3. 
In fact, since the entropy density satisfies s oc —log Q, 
the assumption that Q remains constant is equivalent to 
the hypothesis of an adiabatic inflow, although with a 
different adiabatic constant. 
Using eq. [TU one obtains 
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and replacing the equation above into eq. El one obtains 
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where a Taylor series up to the second order was used. 
The constant k\ in eq. [3] can be calculated by imposing 

0. Thus, one obtains 
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where rioo stands for the dark matter particle number 
density far away from the black hole, where the flow ve- 
locity is negligible. At the critical point, eq. can be 
rewritten as 
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where n c is the number density at the critical point. Re- 
placing eq. [in]into the equation above, one obtains after 
some algebra, 
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which permits to estimate the density at the critical point 
as a function of physical conditions of the halo, far away 
the black hole. Since the flow velocity at the critical point 
satisfies u 2 c « 1/3, from eq. [5] and eq. [THJone obtains to 
the first order 
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and, using the relation above for the density at the crit- 



ical point and the fact that Q Tl 
finally 

M 2 ~ 0.627(^) 



/cr^o, one obtains 



2/3 



(21) 



Replacing this result into eq.JT] one obtains for the critical 
radius 

r c ~ 0.398( — ) 2 '\ g ~ 0.398(— ) 4 / 3 r, (22) 



For typical halos with in the range 50-300 km/s, the 
critical radius is located in the interval 130-40 times the 
horizon radius, well inside the influence sphere. It is in- 
teresting to recall that, for an ideal gas, the critical point 
is located at 
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where is the sound velocity far away from the black 
hole. This relation implies that the critical point is out- 
side the horizon only if 7 < 5/3 (see fl7j|)- 

Once the parameters at the critical point are known, 
the accretion rate is given by 
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where only first order terms were retained. Notice that 
for non-relativistic and non-interacting particles, the ex- 
pected accretion rate is [35j 
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where v<x, is the relative velocity of particles "at infinity" 
with respect to the black hole. This equation predicts 
an accretion rate smaller by a factor 27c 2 /(16v / l25<7 2 ) 
than eq. [M] which is of the order 10 5 , if estimated in the 
core of a typical dark halo. 



III. NUMERICAL APPLICATION 

In this section, an application of the results obtained 
previously is performed. Firstly, we emphasize that we 
do not intend to present here a detailed model for the 



growth of black holes located in the core of halos. This 
problem will be discussed in a future paper. Our purpose 
here is to give only a rough estimate of the dark matter 
fraction accreted by the black hole, assuming that our 
formula (eq. [24]) gives a correct description of the inflow. 

Halos investigated in this work were selected from cat- 
alogs prepared by [33j , which were based on systems iden- 
tified in cosmological N-body simulations, using an adap- 
tive particle-particle/particle- mesh (AP 3 M) code f(3r|). 
The adopted cosmological parameters were h = 0.65, 
= 0.3 and J7a = 0.7, with the power spectrum nor- 
malization taken to be ag = 0.9. Simulations were per- 
formed in a box of side 30/i -1 Mpc, including 256 3 parti- 
cles, corresponding to a mass resolution of 2.05 x 10 s M Q 
and covering the redshift interval < z < 50. Halos were 
initially detected by using a friends-of-friends (FOF) al- 
gorithm ( [37[| ) and, in a second step, unbound particles 
were removed by an iterative procedure. In the total, 
about 9000 gravitationally bound halos were detected at 
z = 0. For more details, see [H, [38|. However, the mass 
history in the redshift range < z < 10 could be followed 
only for 35 halos with enough particles (N > 50). The 
phase space density indicator Q, was calculated accord- 
ing to the prescription detailed in [33| , inside a spherical 
volume whose radius is equal to one tenth of the gravita- 
tional radius of the halo, which is defined as the "core" 
region. 

The present observations indicate that in all massive 
galaxies where the masses of bulges and central black 
holes have been estimated from kinematic data, there is 
a robust correlation between them [H, H3, 53 • More- 
over, other investigations have shown conclusively that 
not only there is a link between the black hole mass and 
the bulge mass but also with the circular velocity [42j . If 
circular motions probe the halo potential, then a relation 
between the present black hole mass and the present halo 
mass is also expected. These relations suggest that black 
holes grew coevally with their host galaxies and associ- 
ated halos [43(. Here, this assumption will be adopted 
and, consequently, the black hole growth rate is given by 
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where a=1.27 is a constant (see, eq.[28lbelow) and Y (t) — 
M b h(t) / M b h(0) is the black hole mass, measured in units 
of its present mass Mbh(0). The right member of eq. [26] 
represents the mass growth rate of the host dark halo. 

At any instant, the total mass of the black hole is the 
sum of the amount of accreted baryonic and dark matter 
masses. If X(t) is fraction of accreted dark matter at 
instant t then, using eq. and the previous definitions, 
one obtains 



dX(t) _ X(t) dY{t) 
dt ~ Y(t) dt 
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where, as in eq. [55] the time is in units of H Q 1 and the 
constant K is equal to K = 27TrG 2 M bh (0)/y / T2EH = 
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FIG. 1: Evolution of the velocity dispersion (km/s) inside the 
core of halos 3 and 211 (upper panel). Evolution of the phase 
density (MQpc~ 3 km~ 3 s 3 ) in core for the same halos (lower 
panel) . Time is given in units of 1 . 

2.143(M Wl (0)/M Q ). The numerical value corre- 
sponds to the phase space density taken in units of 
M Q pc~ 3 km~ 3 s 3 . 

The calculation procedure was the following: for a 
given halo having presently a mass Mjj, the black hole 
mass was estimated by the relation derived by 
namely 

M bh (0) = 1.1 x lO^^^-) 1 ' 27 M (28) 

and, consequently, the value of the constant K can be 
calculated. The system of differential eqs. I2T51 and I2T1 was 
then solved numerically, using the mass and phase space 
density derived from cosmological simulations. The ini- 
tial values of Y and X were varied until the solution gives 
Y = 1 for the present black hole mass. 

For studied halos, with present masses in the range 
7.5 x 10 11 - 9.54 x 10 12 M Q , we have found that the 
fraction of dark matter accreted by the halo is in the 
range 1-10%. 

In order to exemplify, in figure 1 it is shown the evo- 
lution of two halos whose final masses are respectively 



7.46 x 10 n M Q (halo 211) and 8.06 x 1O 12 M (halo 3). 
The upper panel shows the evolution of the velocity dis- 
persion whereas the lower panel shows the evolution of 
the phase space density indicator Q. Both quantities 
were calculated in the core region, as defined above. No- 
tice that to make easier the comparison with observa- 
tions, we have plotted the ID-velocity dispersion and 
that the quantity Q was also calculated with respect to 
uiD- The 3D values can be obtained multiplying the 
ID-velocity dispersion by V3 and dividing Qm by a/27- 
Notice also that the more massive halo (halo 3) has a 
higher velocity dispersion but a lower phase density Q, 
in agreement with the fact that the late quantity, for a 
given redshift, scales as Q oc M -0 82 ([33(]). The black 
hole mass evolution in the considered halos is shown in 
figure 2 (upper panel), where it is possible to see that 
the growth is not uniform, retracing the mass history of 
halos. Sudden increases in mass growth represent merger 
episodes. The fraction of accreted dark matter is shown 
in the lower panel of figure 2. At z ~ 10, the fraction of 
dark matter is about 0.005 for both halos whereas it at- 
tains a value around 0.02 for the less massive halo (halo 
211) and about 0.05 for halo 3. Notice the fast increase in 
the dark matter contribution for z < 2, consequence of a 
decreasing merger rate after such a redshift. This affects 
the dark matter accretion in two aspects: firstly, decreas- 
ing the contribution of baryonic matter to the black hole 
growth and secondly, the decreasing merger rate implies 
also that Q decreases more slowly than the black hole 
growth rate and, consequently, the relative contribution 
of dark matter increases. 

In figure 3 the baryonic and dark matter accretion 
rates are compared for both halos. Different peaks in 
the baryonic accretion rate are associated to minor or 
major merger episodes, with the maximum rate occur- 
ring around z ~ 1.5. Halo 3 is an example where for 
two short time intervals the accretion rate of dark mat- 
ter became higher than the baryonic accretion rate. The 
baryonic accretion rate will be probably reduced by a 
factor of 2-6 if the star formation process is included, an 
aspect presently under investigation, based on detailed 
cosmological simulations. 



IV. CONCLUSIONS 

The accretion of dark matter into black holes located 
in the center of galaxies and their halos was studied. The 
accretion rate was derived by solving the equations for a 
relativistic steady spherical symmetric inflow, under the 
assumption that the phase space density indicator Q = 
pjo 3 remains constant. This hypothesis is equivalent to 
say that the inflow is adiabatic, since the entropy density 
satisfies s oc — log Q. 

We found that our solution predicts that the critical 
radius is always outside the black hole horizon and that 
the flow radial velocity at the critical point is of the order 
of few percent of the speed of light. After the crossing the 
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FIG. 2: In the upper panel it is shown the mass evolution 
of black holes while in the lower panel the evolution of the 
fraction of accreted dark matter is shown. 

critical point, numerical solutions of equations describing 
the flow indicate that the velocity and the density have 
power-law profiles, e.g., u oc 1/r 6 and p oc 1/r 1 ' 4 , a 
behavior quite close to that resulting from the classical 
"free-fall" solution. Notice that the growth of black holes 
embedded in a collisionless fluid was also studied by con- 
sidering that orbits shrink inside the influence sphere but 
conserving their adiabatic invariants [3] ■ In this case, a 
density "spike" with a profile p oc 1/r 1 - 5 results, similar 
again to our inflow model. 

The derived accretion rate for a steady flow is propor- 
tional to the square of the black hole mass and to the 
dark matter phase space density indicator Q = poo/a^. 
For black holes hosted by typical dark halos, the present 



accretion rate is about five orders of magnitude higher 
than the usual rate derived by using the capture cross 
section for non-relativistic and non-interacting particles. 

Numerical applications of the derived formula were 
made by computing the growth of black holes inside ha- 
los issued from cosmological simulations. We found that 
dark matter contributes to no more than 10% of the to- 
tal mass accreted by black hole seeds. This result is con- 
sistent with the fact that the time integrated bolometric 
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FIG. 3: Evolution of the baryonic and dark matter accretion 
rates for halo 211 (upper panel) and halo 3 (lower panel). 



luminosity of quasars gives a present black hole mass den- 
sity comparable with other independent estimates, indi- 
cating that most of the accreted mass by seeds is baryonic 
in origin. 
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